Residence time statistics for iV blinking quantum dots and other stochastic processes 
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We present a study of residence time statistics for TV blinking quantum dots. With numerical 
simulations and exact calculations we show sharp transitions for a critical number of dots. In 
contrast to expectation the fluctuations in the limit of TV — » oo are non-trivial. Besides quantum 
dots our work describes residence time statistics in several other many particle systems for example 
TV Brownian particles. Our work provides a natural framework to detect non-ergodic kinetics from 
measurements of many blinking chromophores, without the need to reach the single molecule limit. 
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Residence time statistics has attracted considerable in- 
terest as a tool describing non-equilibrium phenomena. 
For a fluctuating time trace I(t) the residence time t + 
is the total time the signal remains beyond some thresh- 
old. A classical result in this field is Levy's arcsine law 
0, 0. Then t+ is the residence time of a one dimen- 
sional Brownian particle in half space. In contrast with 
naive expectation the PDF of t + /t (t is the measure- 
ment time) is not centered on its mean, instead it has a 
U shape with singularities described by the persistence 
exponent a = 1/2 for simple diffusion [3| (see details be- 
low). Similar bi-modal distributions were obtained for 
random walks in disordered environments NL the persis- 
tence problem of diffusion in dimension d |5j, a melting 
heteropolymer Sinai model @ and for the fluorescence 
signal from a quantum dot (QD) [7, 8] . A single QD when 
interacting with a continuous wave laser field exhibits 
blinking with power law kinetics which is related to non 
trivial residence times, non-stationary and non ergodic 
behavior [7H1C|. At random times the dot will switch 
from a bright state where many photons are emitted to 
a dark state (see Fig. [T]). Probability density function 
(PDF) of on and off times are power laws ip{r) ~ r~( 1+Q ) 
at least for low laser intensity [ll| . In most cases a ~ 1/2 
though 1/2 < a < 1 was also reported. Most interest- 
ingly the average on and off times (t) = J °° tiP{t)At 
diverges, i.e. the dynamics is scale free. Several physical 
models explaining this behavior were sugg ested [9|, llfj| . 
for example based on normal diffusion 8] , however so far 
no experimental smoking gun provides confirmation to a 
specific mechanism responsible for the power law blink- 
ing. Further the effect is not limited to QDs, it is found 
also for organic single molecules provided that they are 
embedded in a disordered material [Io| . One can rightly 
wonder why after many decades of spectroscopy we find 
surprises once individual objects are detected [12J. Can 
we not infer the strange blinking from measurements of 
many dots/single molecules? In particular can the bro- 
ken ergodicity previously measured for an individual dot 
@, H| be detected also for an ensemble of dots? 

For that aim we investigate the problem of residence 
time statistics for TV blinking dots. We will soon show 



FIG. 1: Model two state blinking QD dot with a = 1/2. For 
TV = 1 (lower panel) power law distributed sojourn times in 
bright +1 and dark states lead to long residence times, of 
the order of the measurement time, in both states. Upper 
panel exhibits a superposition of the signal from TV = 501 
processes. In this Letter we investigate the fluctuations of 
the residence time, namely the time the signal remains above 
the ensemble mean (the dashed line) showing they are not 
trivial even in the limit of TV — > oo. 



that this problem can be mapped onto an interesting 
many particle Brownian problem and other stochastic 
processes. Surprisingly non trivial residence time statis- 
tics is found here even in the limit of TV — > oo . This is an 
important indication that by the analysis of light emit- 
ted from an ensemble of chromophores we may detect 
the strange kinetics without the need to reach the single 
molecule limit. We investigate transitions in the PDF 
of residence times (from U to n shape, as explained be- 
low) which are controlled by the number of dots TV. We 
show that a critical number of dots marks a border line 
between the single particle domain and the many parti- 
cle limit. Surprisingly we find values of a below which 
U shape PDFs are robust, namely for certain a the U 
shape PDF is found even for TV — > oo. Roughly speaking 
this is an important indication that elements of Levy's 
arcsine law survive the TV — > oo limit, hence persistence 
exponents [3| from signals composed of a super-position 
of many elements are more common than expected. Our 
main analytical calculation is an exact expression for the 
variance of the residence time which is used to identify 
transitions between sub and super uniform statistics de- 
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fined below. 

Model and observable. We consider A independent 
though statistically identical QDs. Each dot undergoes 
the following simple renewal process. The dot can be in 
a bright state where the intensity of light is Ii (t) = 1 or 
a dark state with intensity Ii(t) = 0. Sojourn times in 
state on and off are independent identically distributed 
random variables with a common PDF V'( r )- V'( T ) is 
moment less, namely for large r ip(r) ~ Ar~( 1+a ^ and 
< a < 1. The signal of A dots is l N {t) = Y^Li M*) 
whose ensemble mean is obviously (Xjv) = A/2. Here 
we investigate the time t+ the process remains above its 
mean (see Fig. [J). The residence fraction is defined as 
< p + = t + /t < 1 where t is the total measurement 
time. 

Relation with Brownian motion. The model and its 
variants describe many physical systems and processes 
beyond QDs e.g. well known Levy walks which de- 
scribe tracer diffusion in turbulent flow, certain chaotic 
systems [3] and recently models of 1// noise [lil [lil]. 
Probably the best well known example is Brownian mo- 
tion. Consider a single Brownian particle in one dimen- 
sion and let Ii(t) = 1 for a particle being in x > and 
h(t) = for x < 0. As well known the PDF of first 
passage times and hence of sojourn times in x < or 
x > follow power law behavior ip( T ) ^ t -3 / 2 and hence 
a = 1/2. For a single particle the time t + is simply the 
residence time in half space x > which follows the clas- 
sical arcsine distribution of P. Levy (sec details below). 
For TV Brownian particles t + is the time the majority of 
A particles are on x > 0. Naively one would expect that 
this time in the A — > oo limit would be equal to half of 
the measurement time. This turns out to be wrong and 
features of Levy's arcsine law survive the N — > oo limit. 

Known results for a single QD A — 1. Lamperti ob- 
tained the PDF of the residence fraction < p + < 1 



la( P + ) = 



[sin (ira) /tt] (p + ) a - 1 {l - p + ) c 



(p + ) 2a + 2(p+) a (l - p+) a cos(Tra) + (1 - p+) 2a ' 

(1) 

When a = 1/2 the mentioned arcsine PDF is obtained 
which has a U shape. When a — > 1 the PDF Eq. ^ 
approaches a delta function on its mean (p + ) =1/2 which 
corresponds to ergodic behavior. In the opposite limit 
a — > we have two delta functions on p + = and p + = 1. 
Then the QD is either in the dark state or the bright 
state for the whole duration of the measurement. An 
important feature of the Lamperti PDF is that for < 
a < 1 the maximum of the PDF is on these rare event: 
the PDF diverges on p + = and p + = 1. This clearly 
reflects sojourn times in either the on or the off state 
which are of the order of the measurement time. 

The PDF of the residence fraction was obtained numer- 
ically and is presented in Fig. [5]for a = 0.57 and a = 0.5. 
For a = 0.57 a transition from a U shape PDF for < 5 
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FIG. 2: Left panel a — 0.57 the PDF of residence fraction 
exhibits a transition between a U to n shape PDF which is 
controlled by N . In contrast such a transition is not observed 
for a — 0.5 and for smaller values of a (not shown). 



to a n shape PDF for N > 5 is found, while N — 5 gives 
a uniform PDF. In contrast, for a — 0.5 the PDF for 
all N has a U shape. These findings indicate: (i) even 
in the large N limit the fluctuations of residence times 
remain non-trivial, namely we do not see a convergence 
of the PDF to a delta function and (ii) that there exists a 
critical a below which the PDF of residence fractions has 
always (for any N) a U shape. However, simulations are 
limited to finite N and finite measurement time, hence 
we now turn to analytical theory. 
The residence fraction for N QDs is 



P 



[l N (t')-N/2] df 



(2) 



where 6{x) is the step function: 9(x) — 1 if its argument 
is positive, otherwise it is zero. For convenience we will 
consider only odd N. The states of the system are de- 
termined by a vector \Ii(t), ...Ii(t)...I N (t)]. States con- 
tributing to the considered residence fraction are those 
with at least (N + 1)/2 processes in bright state ij(i) = 1 
and we have 2 N /2 such states. We can rewrite the resi- 
dence fraction as a sum over time averages of these states 
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(3) 



Here H^ =1 Ii»^ is a sum over all permutation of the prod- 
uct Ii...In with j terms of the type 1 — Ii and and 
N — j terms of the type Ii. For example for A = 3 
we have Ilf =1 Ii» = hhh for j = and H^^. 1 = 
(1 - h ) h h + h (1 - h ) h + h h (1 - h ) ■ This means that 
p + has contributions from states (1, 1, 1) (all dots bright) 
or (0, 1, 1), (1, 0, 1), (1, 1, 0) since for these states the total 
intensity Z}v=3(t) is above its mean (X)n=3 = 3/2. 

Sub and super uniform statistics. It is easy to see that 
the mean is (p + ) = 1/2 which is expected from symme- 
try. More interesting is the variance a 2 = ((p + ) ) — {p + ) 2 . 
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TABLE I: Variance of the residence fraction a 2 for various 
AT. 



The variance can be used to quantify different types of 
fluctuations. If a 2 > 1/12 we define the fluctuations as 
super- uniform while a 2 < 1/12 is sub-uniform. Distribu- 
tions of residence fractions with U (n) shape are super- 
uniform (sub-uniform) respectively. Squaring Eq. (|3]) 
and averaging with respect to the random process 
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±J2a n (i) fdtx [* dhintjii - i(t 2 )]) i {i( tl )i(t 2 )) 
1=0 J ° J ° 

(4) 

We see that the fluctuation of the residence frac- 
tion is determined by the intensity correlation function 
(I(h)I(t 2 )) (here (1(h) [1 - I(t 2 )]) = 1/2- (J(ti)/(t a )». 
A lengthy combinatorial calculation, which we will pub- 
lish elsewhere, yields the coefficients An (I) 



N- 



A N (l) 



(^)2^ 1 [l-S N (l)] if/<^ 



(5) 



(7)2^-^(^-0 ifi> 



N-l 
2 ' 



where 17 



S N (l) = 2 £ 



i=0 



l\ 1 1 



i 2 l 2 l - 



E 



N - A 1 1 

i J 2^— I— i 2* 



(6) 

This representation will soon be useful since we identify 
binomial expansions in it. 

Correlation functions. Many natural processes exhibit 
stationary behavior namely (I(t\)I(t 2 )) is a function of 
the time difference \t 2 — t\ \ . However many other systems 
exhibit aging 



(I(t 1 )I(t 2 ))=C(t 1 /t 2 ) h/t 2 <l, 



(7) 



which is common for a system with infinite mean waiting 
time. For the two state QD process [l], H[ 



C(z) = - 



sin 7tcv 



B(z; a, 1 — a) + 1 



< z < 1 (8) 



and B(z] a, 1 — a) — J Q x a 1 (l — x) "dec is the tabulated 
incomplete Beta function. 
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FIG. 3: The variance of the residence fraction versus a for 
various number of processes N. When a — >• 1 the process 
turns ergodic and the variance is zero. For a < 1 even when 
iV — i> 00 the fluctuations remain finite, indicating that in prin- 
ciple the non-ergodic kinetics so far detected on the single 
particle level can be found already for a large ensemble of 
. particles. Even more surprisingly, for N — s- 00 the fluctua- 
tions remain super-uniform, i.e. the variance is larger than 
1/12, when a < a c — 0.518.... This implies that a n shape 
PDF of the residence times does not describe fluctuations 
when a < a c . 



Fluctuations for finite N. Inserting Eq. ([SJ in Eq. Q 
changing variables to z = ti jt 2 we find exact expressions 
for the fluctuations. For N — 1 a 2 = 1 — a, while for 

N = 3 



o~ 2 = 4C 3 



3C 2 



V 1 



where 



C n = I C n ( Z )dz. 
'0 



(9) 



(10) 



Cases N = 5, 7, 9 are reported in Table 1 and similarly we 
obtain a 2 for any finite N using the exact expression for 
An (I). With Table 1 we can determine whether fluctua- 
tions are sub or super uniform. For example for N = 5 
we find uniform statistics a 2 = 1/12 when a ~ 0.57. 
This perfectly matches the simulations in Fig. [5] which 
exhibit a uniform distribution of the occupation fraction 
for these parameters. For a = 1/2, i.e. Brownian motion, 
Eqs. (j8TTO| give 



m=0 



(nr)- 



r (1 + m, — in) — r (1 + rn, ii 
(11) 



For a ^ 1/2 we solve the integral Eq. (ITUl) with Math- 
cmatica [18j and with Table 1 obtain numerically exact 
expressions for a 2 . In Fig. [3] we show a 2 versus a for 
various AT. When a = 1 we have C n = (I) 2n = 4~ n 
and then using Table 1 we get a 2 = which corresponds 
to the ergodic phase. In the opposite limit a — > we 
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have C n = (I) n = 2~ n and then a 2 = 1/4. This is 
the expected behavior since when a-}0a single parti- 
cle remains in one state (either bright or dark) for the 
whole measurement time and hence with probability 1/2 
the signal composed of an ensemble of an odd number of 
particles is above its mean all along the measurement. As 
shown in Fig. [3] the fluctuations are the strongest when 
N = 1 which is expected. Surprisingly, we see that as the 
number of particles is increased the fluctuations remain 
finite. This means that the non ergodic behavior found 
for a single QD 0, [|| can be detected also for a large 
ensemble of dots. A measurement of the residence time 
t + remains a random variable even when N is large and 
a single measurement is not a reproducible result. We 
now turn to find the variance a 2 in the N — > oo limit. 

Limit N — > oo. The limit of large N and I is taken 
in such a way that x = l/N remain finite. We then use 
Stirling's formula to approximate the double binomial 
sum in Eq. ([6]) with a double integral over a pair of 
Gaussians. The first integration yields an error function 
and we find 



S N (l) 




xN _z«£^ , / v^; 
e 2 Erf _ _ 



Nxz 



dz. (12) 



From a table of integrals Sn(1) ~ ^ArcSin(y / x) and 
hence we find 



A N (I) 



>JV-1 



1 ArcSin 

7T 




(13) 



super-uniform statistics and surprisingly this property is 
maintained also for a large ensemble of particles. Insert- 
ing Eq. © in Eq. ([T5J we find that if a < a c = 0.518 • • • 
the fluctuations are super-uniform in the limit N — > oo 
(see Fig. This theoretical prediction is in agreement 
our simulations presented in Fig. [5] which shows super- 
uniform (sub-uniform) behaviors when a < a c (a > a c ) 
respectively. 

Bouchaud [2(J, introduced the profound concept of 
weak ergodicity breaking in the context of glassy dynam- 
ics. It holds for systems with power law sojourn times 
and with a phase space not broken into mutually inacces- 
sible regions, e.g. blinking QDs. For a macroscopic sys- 
tem like a glass Bouchuad assumed that measurements 
are made on a large number of sub-systems (micro spin- 
glasses of clusters of magnetic atoms in his case) and that 
such measurements are disorder averaged. He claimed 
that performing dynamical experiments on mesoscopic 
samples, where N is small, would yield irreproducible 
results. Indeed we thoroughly believed this scenario at 
the start of this project, namely our naive expectation 
was that as the number of QDs increases we will ap- 
proach an ergodic behavior. However, it turns out that a 
measurement of the residence time of N — > oo QDs is ir- 
reproducible. Weak ergodicity breaking can be detected 
for a large ensemble of sub-systems. 
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Inserting this expression in Eq. ^) and using Eq. (|7J| we 
find (after change of variables to z = t-y/t-i) 
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; ArcSin J- ( , ) [1 - 2C(z)f [2C(z)\ 
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In the large N limit the binomial part of the integral: 

is narrowly centered around 
its mean [1 — 2C(z)] N. Hence we find 



(7) [l-2C(z)] c [2C(z)] 



N-l 



lim a 

N— >co 



1 



dzArcSinv/2C(z) - -, 



(15) 



which is the main equation of this manuscript. 

Eq. (fT5j) shows that even when N — > oo the statistics 
of blinking QDs is non ergodic: the fluctuation of the 
residence time arc finite even for a large ensemble. Only 
in the limit a — > 1 we do get expected ergodic behavior 
(T 2 = 0. In Fig. [3] we plot Eq. (fT"5|) versus a and classify 
sub-uniform and super-uniform statistics. For the Brow- 
nian case a — 1/2 we get limjv_ i . 00 a 2 = 0.086 • • • which 
is super-uniform, hence this rigorously shows that a n 
shape PDF is not found Levy's arcsine law exhibits 
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